The Atacama Cosmology Telescope: Cross-Correlation of CMB Lensing and Quasars 



(N 

o 

OV 



O 

u 

6 

o3 



> 
o 

(N 



13 



Blake D. Sherwin, 1 ^ Sudeep Das, 2,1 ' 3 Amir Hajian, 4 Graeme Addison, 5 J. Richard Bond, 4 Devin 
Crichton, 6 Mark J. Devlin, 7 Joanna Dunkley, 5 Megan B. Gralla, 6 Mark Halpern, 8 J. Colin Hill, 3 
Adam D. Hincks, 4 John P. Hughes, 9 Kevin Huffenberger, 10 Renee Hlozek, 3 Arthur Kosowsky, 11 

Thibaut Louis, 5 Tobias A. Marriage, 6 ' 3 Danica Marsden, 12 Felipe Menanteau, 9 Kavilan Moodley, 13 
Michael D. Niemack, 14, 1 Lyman A. Page, 1 Erik D. Reese, 7 Neelima Sehgal, 3 Jon Sievers, 1 
Cristobal Sifon, 15 David N. Spergel, 3 Suzanne T. Staggs, 1 Eric R. Switzer, 16 ' 1 and Ed Wollack 17 

1 Dept. of Physics, Princeton University, Princeton, NJ, USA 08544 
2 BCCP, Dept. of Physics, University of California, Berkeley, CA, USA 94720 
3 Dept. of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ USA 08544 
4 CITA, University of Toronto, Toronto, ON, Canada M5S 3H8 
° Dept. of Astrophysics, Oxford University, Oxford, UK 0X1 3RH 
6 Dept. of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218-2686 
7 Dept. of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA, USA 19104 
8 Dept. of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada V6T 1Z4 
3 Dept. of Physics and Astronomy, Rutgers, The State University of New Jersey, Piscataway, NJ USA 08854-8019 
10 Department of Physics, University of Miami, Coral Gables, Florida 33146 
11 Dept. of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA, USA 15260 
12 Department of Physics, University of California, Santa Barbara, CA 93106 US 
13 Astrophysics and Cosmology Research Unit, Univ. of KwaZulu-Natal, Durban, 404^ South Africa 
14 NIST Quantum Devices Group, 325 Broadway Mailcode 817.03, Boulder, CO, USA 80305 
15 Departamento de Astronomia y Astrofisica, Facultad de Fisica, 
Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 
16 Kavli Institute for Cosmological Physics, 5620 South Ellis Ave., Chicago, IL, USA 60637 
17 Code 553/665, NASA/Goddard Space Flight Center, Greenbelt, MD, USA 20771 

We measure the cross-correlation of Atacama Cosmology Telescope CMB lensing convergence 
maps with quasar maps made from the Sloan Digital Sky Survey DR8 SDSS-XDQSO photometric 
catalog. The CMB lensing-quasar cross-power spectrum is detected for the first time at a significance 
of 3.8a, which directly confirms that the quasar distribution traces the mass distribution at high 
redshifts z > 1. Our detection passes a number of null tests and systematic checks. Using this 
cross-power spectrum, we measure the amplitude of the linear quasar bias assuming a template for 
its redshift dependence, and find the amplitude to be consistent with an earlier measurement from 
clustering; at redshift z ~ 1.4, the peak of the distribution of quasars in our maps, our measurement 
corresponds to a bias of b = 2.5 ± 0.6. With the signal-to-noise ratio on CMB lensing measurements 
likely to improve by an order of magnitude over the next few years, our results demonstrate the 
potential of CMB lensing cross-correlations to probe astrophysics at high redshifts. 



I. INTRODUCTION 

As the photons of the cosmic microwave background 
(CMB) travel through the universe, they are gravita- 
tionally deflected by the web of matter through which 
they pass. In the CMB sky we observe today, these 
deflections are imprinted as arcminute-scale distortions 
of small scale temperature fluctuations [l|, [2| . The mi- 
crowave background thus contains not only information 
about the primordial universe at redshift w 1100, but 
also about the matter density fluctuations in the lower 
redshift universe. 

The lensing deflection held d(n) points from the direc- 
tion n in which a CMB photon was received to the direc- 
tion from which it was emitted. This deflection field can 
be determined from the measured CMB because lensing 
changes the statistics of small scale unlensed CMB fluc- 



" bsherwin@princeton.edu 



tuations in a characteristic way, introducing correlations 
between different Fourier modes. By measuring corre- 
lations between pairs of Fourier modes that would be 
uncorrclated in the absence of lensing, one can estimate 
d [3( and hence the lensing convergence k = — V-d/2 
(a useful quantity because it is a direct measure of the 
projected matter density, see Eq. 1). Using a quadratic 
estimator, lensing was first measured in cross-correlation 
with radio sources and galaxies by 0, [H (using WMAP 
data) and in auto-correlation by 6] (using Atacama Cos- 
mology Telescope data). More accurate measurements of 
both the lensing power spectrum and the lensing-galaxy 
cross-correlation were recently reported by the South 
Pole Telescope 

CMB len sing measurements are a powerful cosmolog- 
ical probe [1C§ because they are sensitive to both the 
growth of density fluctuations and the geometry and size 
of the universe, yet are relatively insensitive to both 
instrumental and astrophysical systematic errors [f| 0]- 
Lensing measurements can constrain the properties of 
dark energy 18|, the amplitude of density fluctuations 
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and the masses of neutrinos [11] . They can also con- 
strain the properties of biased tracers of the matter dis- 
tribution. The focus of this paper is the cross-correlation 
of CMB lensing maps with one such tracer - quasars. 

Quasars are among the most luminous objects in the 
universe. Their immense luminosity is believed to be 
powered by accreting supermassive black holes [HI, [l3| 
which reside at the center of almost every massive galaxy 
[l~4| . As the activity of quasars and the rate of star for- 
mation appear to be linked [l||, they are a crucial ele- 
ment in our present understanding of galaxy evolution. 
Measurements of the relation between dark matter and 
the distribution of quasars can inform us about quasar 
properties such as the masses of the dark matter halos 
that host the quasars, the scatter in the quasar-halo mass 
relation and the quasar duty cycle (see e.g. [16|). Such 
measurements of quasar properties will, in turn, improve 
our understanding of structure formation and galaxy evo- 
lution. 

Both the number density of quasars and the strength 
of CMB lensing in a certain direction depend on the pro- 
jected dark matter density in this direction, and quasars 
are most common at the redshifts that produce the 
largest lensing deflections. This implies that the CMB 
lensing and quasar fields should be strongly correlated 
17] . Measuring the cross-power spectrum and comparing 
it to theoretical calculations, we can determine the pro- 
portionality factor which relates a fluctuation in matter 
density to a fluctuation in quasar number density. This 
proportionality factor is known as the quasar linear bias 
b (defined as b = 5 q /S where S q and 8 are the fractional 
spatial overdensities of quasars and matter respectively). 

In this work we present the first measurement of the 
CMB lcnsing-quasar cross-power spectrum, and use it to 
derive a constraint on the quasar bias. The paper is or- 
ganized as follows. Section II presents the theoretical 
background underlying the CMB lensing-quasar cross- 
correlation. Section III explains how the lensing and 
quasar maps used in our analysis are constructed, and de- 
scribes the resulting cross-spectrum measurement. The 
constraint on quasar linear bias we obtain from the cross- 
power spectrum is presented in section IV. Null tests and 
systematic checks are discussed in section V. All calcula- 
tions assume a ACDM cosmology with WMAP5 param- 
eters pj] and cr 8 = 0.819. 



II. THEORETICAL BACKGROUND 



The relevant kernel (assuming a flat universe) is 
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where r/(z) is the comoving distance to redshift z, n is 
a direction on the sky, rj LS is the comoving distance to 
the last scattering surface, zls is the redshift of the last 
scattering surface, H{z) is the Hubble parameter, and fio 
and Hq represent the present values of the matter density 
parameter and the Hubble parameter, respectively. 

The fractional overdensity of quasars in a direction n 
is given by q(n), which - assuming a linear bias relation 
between the distribution of matter and quasars - is given 
by 



q(n)= dzW q (z)S(r)(z)h,z), (3) 
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where the kernel is 
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(see pjj) ■ Here b is the linear bias, dN/dz is the red- 
shift distribution of quasars, and s is the variation of 
the number counts of quasars N(< m) with the limit- 
ing magnitude m at the faint limit of the quasar catalog, 
s = dlog w N/dm. The second term in Eq. 4 is the magni- 
fication bias: the change in source density in a certain di- 
rection induced by lensing magnification. For the quasar 
catalog used in this work (which has s « 0.043) this term 
is significantly smaller than the first term (w 15% of its 
magnitude) and is negative. 

Combining Eqs. 1 and 3 gives the expected lensing- 
quasar cross-power spectrum in the Limber approxima- 
tion [l^: 

CT = J ^^W K (z)W\z)P(k = £/ V (z),z), (6) 

where P(k, z) is the matter power spectrum at wavenum- 
ber k and redshift z. 



Cosmological weak lensing effects can be described us- 
ing the convergence field k, which is equal to a weighted 
projection of the matter overdensity 5 [l| 

k(&)= dzW K (z)S(r](z)h,z). (1) 

Jo 



III. CROSS-CORRELATING CMB LENSING 
AND QUASARS 

A. The ACT CMB Lensing Convergence Maps 

The lensing convergence fields used in our analysis are 
reconstructed from CMB temperature map s made by the 
Atacama Cosmology Telescope (ACT) [23M25j | , a 6m tele- 
scope operating in the Atacama desert in Chile. These 
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CMB temperature maps are obtained from observations 
made during 2008-2010 in the 148 GHz frequency band 
and calibrated as in 26] . The maps consist of six patches, 
each of size 3x18 degrees, in the Sloan Digital Sky Survey 
(SDSS) Stripe 82 region [23], with a map-averaged white 
noise level of 21 /iK arcmin. Radio and IR point sources 
as well as Sunyaev-Zel'dovich (SZ) clusters detected with 
a matched filter at a signal-to-noise ratio greater than 5 
as in [6J are masked and inpainted using the methods of 
1- 

The convergence maps are reconstructed from the 
CMB temperature maps as in [f| using the minimum vari- 
ance quadratic estimator procedure described in [3j . The 
estimator works as follows. While the unlensed CMB is 
statistically isotropic, any lensing deflection imprints a 
preferred direction into the statistical properties of the 
CMB. This corresponds mathematically to the fact that 
formerly uncorrelated modes of the isotropic unlensed 
CMB temperature field are correlated by lensing, with 
the correlation proportional to the lensing deflection. We 
can hence estimate the lensing convergence by measur- 
ing the correlation between pairs of Fourier modes using 
a quadratic estimator: 



k(L)=N(L) / d 2 l/(L,l) T(1)T(L-1). 



(7) 



Here 1, L are Fourier space coordinates, N is the nor- 
malization function (which ensures that the estimator 
is unbiased) and / is a weighting such that the signal- 
to-noise ratio on the reconstructed convergence is maxi- 
mized (see [3J] for the details of these functions). The final 
normalization is obtained from Monte-Carlo simulations 
(applied as a small correction to a first-approximation 
normalization which uses the data power spectrum); the 
simulated maps used in this Monte- Carlo calculation are 
constructed to match ACT data in both signal and noise 
properties as in [6]. In our lensing reconstruction, we 
use a wider range of scales in the temperature map than 
in Q, filtering out modes below I = 500 and above 
I = 4000. As can be seen in Eq. 7, lensing informa- 
tion at a scale £ is obtained from two temperature modes 
separated by £, so that this filtering does not preclude us 
from obtaining small-^ lensing modes. Intuitively, this is 
because we deduce the distribution of large scale lenses 
from the distortion of small scale temperature fluctua- 
tions. 

We subtract from the reconstructed ACT convergence 
maps a simulated mean field map (k), obtained from 
480 realizations of simulated reconstructed lensing maps. 
This map is non-zero due to correlations induced by win- 
dow functions and noise which, to the convergence esti- 
mator, appears as a small spurious lensing signal which 
must be subtracted. 




FIG. 1. The redshift distribution of SDSS quasars used to 
construct our maps of fractional quasar overdensity, normal- 
ized to a unit maximum. The corresponding redshift bins are 
shown with blue filled circles; they are interpolated to give the 
continuous curve used in our theory calculations (blue dashed 
line). For comparison, the red dotted line shows the lensing 
kernel W K (z), again normalized to a unit maximum. 



B. The SDSS Quasar Maps 

In this work we use the SDSS-XDQSO photometric 
quasar catalog [111, [22j , extracted from the SDSS Data 
Release 8. The analysis used in this catalog separates 
the population of quasars and foreground stars using a 
probabilistic model in flux space. This analysis assigns a 
probability of being a quasar to every point-source with 
de-reddened i-band magnitude between 17.75 and 22.45 
mag in the SDSS imaging. Though the catalog extends 
out to z > 4, we do not use the highest-redshift sources 
with redshifts z > 3.5. This reduces the shot noise error 
on our measurement, as there are very few such sources. 
Using this catalog, we construct a map of the fractional 
overdensity q of quasars across the same 324 square de- 
grees on which we perform our lensing reconstruction. 
We include in our quasar maps, with unit weight, all 
point sources with a greater than p = 0.68 probability of 
being a quasar. As the probability distribution is non- 
uniform, the residual level of stellar contamination can 
be calculated from the catalog probabilities to be 6% (we 
discuss later how this is accounted for in our theoretical 
calculations). A spectroscopic quasar sample at high red- 
shifts (2.2 < z < 3.5) does show a contamination fraction 
which is « 15% larger than that predicted by the cata- 
log (see [Hj]); however, the majority of the quasars we 
consider are at lower redshift where the XDQSO model 
estimates should be significantly more accurate (high red- 
shift quasar selection is less accurate because at z s» 2.8 
the quasar and stellar loci cross in color space). We 
neglect the error on the stellar contamination fraction, 
which is any case smaller than the statistical error (and 
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would only reduce the measured signal, not increase it). 
The area covered by our quasar maps contains on average 
75 quasars per square degree. The redshift distribution 
of the quasars in our maps is shown in Fig. 1, along with 
the CMB lensing kernel, which has a very similar redshift 
distribution. 



C. The CMB Lensing - Quasar Cross- Power 
Spectrum 

The cross-power spectrum of the ACT CMB lensing 
maps and the SDSS quasar maps is shown in Fig. 2. 
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FIG. 2. The CMB lensing - quasar density cross-power spec- 
trum, with the data points shown in blue (the covariance be- 
tween different data points is negligible). The significance of 
the detection of the cross-spectrum is 3.8cr. The green solid 
line is a theory line calculated assuming the fiducial bias am- 
plitude. This theory line is reduced by 6% to account for the 
expected level of stellar contamination of the quasar sample. 

The error bars on the data points are calculated theo- 
retically as proportional to -J C^Cf 1 , with an additional 
factor calculated from the number of independent pixels 
in the bin corresponding to each data point. The spectra 
used in this calculation are the full data spectra which 
include both signal and noise (including Poisson noise). 
Bootstrap error estimates from splits of our data are con- 
sistent with this calculation. For comparison, we also 
calculate error bars using Monte-Carlo methods, cross 
correlating 480 realizations of simulated reconstructed 
lensing maps (containing signal and realistic noise) with 
the quasar data maps. Both methods assume the two 
maps are uncorrelated; such calculations are very good 
approximations to the true error on the cross-correlation 
because both maps are noisy so that C^CJ 9 » (C^ q ) 2 . 
The error bars obtained from Monte-Carlo methods are 
nearly indistinguishable from the theory error bars, and 
lead to the same detection significance. (We also verified 
that replacing the quasar data maps with 480 simulated 
maps with the same number of randomly distributed 
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FIG. 3. Blue dashed line: fiducialquasar bias template (inter- 
polated from the data points of [3l[), used in the theoretical 
calculation of the CMB lensing - quasar cross-power spec- 
trum. Red solid line: the bias amplitude (b/b&d = 1.02) best 
fit by the measured cross-power spectrum; red dashed lines: 
the ±1<7 error ranges of this amplitude. Note that only one 
degree of freedom is constrained: the overall bias amplitude 
of an assumed redshift evolution. 

sources leads to similar, though slightly smaller, errors.) 
The Monte-Carlo estimates of the errors also allow us to 
calculate the full covariance matrix. The off-diagonal ele- 
ments are negligible compared to the diagonal elements; 
for every bin, the covariance between neighboring bins 
was less than 4% of the bin autocorrelation. We thus 
neglect covariance between different data points in our 
analysis. 

Also shown in Fig. 2 is a theoretical calculation of the 
expected cross-power spectrum obtained from Eq. 6. In 
this calculation the matter power spectrum was com- 
puted using the CAMB software [29|. The non-linear 
(HALOFIT, 30]) matter power spectrum was used; how- 
ever, using a linear matter power spectrum instead only 
slightly changed the computed cross-spectrum (as most 
of the signal arises from angular scales corresponding to 
linear scales in the matter power spectrum, where the lin- 
ear and non-linear matter power spectra hardly differ). A 
fiducial bias model for this calculation is obtained by in- 
terpolating the central measured bias values of [31] (aver- 
aging the values obtained with and without the inclusion 
of negative points in the correlation function). The frac- 
tional error on these central bias values is below 10% at 
low redshifts z < 2, but rises to ~ 20% at z ~ 4 — 5. The 
values were obtained from measurements of the ampli- 
tude of the quasar correlation function (which is sensitive 
to the bias) for an SDSS spectroscopic quasar sample. 
This fiducial bias model is shown in Fig. 3. As the the- 
oretical cross power spectrum does not depend strongly 
on the detailed form of the bias model, we use this mea- 
surement as a convenient fiducial template, though the 
spectroscopic catalog used in this measurement does not 
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extend to as faint a magnitude as the photometric cata- 
log we use to make quasar maps. Despite this, the quasar 
power spectrum predicted by this fiducial bias model is 
consistent with the power spectrum of our quasar maps. 
The calculated theoretical cross-power spectrum is re- 
duced by 6% to account for stellar contamination; while 
stars are uncorrelated with lensing, they contribute to the 
average density of sources, and so cause us to calculate a 
fractional quasar overdensity that is 6% too small. 

The data fit the theory curve (which assumes the fidu- 
cial bias model) well, with a chi-squared value for this 
curve of Xt heory = 13.2 for 10 degrees of freedom. We ob- 
tain the significance of our detection of the cross-power 

spectrum by calculating yxLn ~ XiLory where Xnull is 
calculated for the null line. The significance of the detec- 
tion is found to be 3.8<r. 

IV. A CONSTRAINT ON THE QUASAR BIAS 

We calculate a constraint on the linear bias of quasars 
from the lensing-quasar cross-power spectrum. To do so 
we assume a bias template of the fiducial shape shown 
in Fig. 3, but rescaled by a constant factor for all red- 
shifts. We calculate the likelihood as a function of this 
scaling parameter 6/6ga and plot it in Fig. 4. Our re- 
sult, 6/6fld = 1-02 ± 0.24, is consistent with the fiducial 
bias model, i.e. a value of unity. (Due to the small neg- 
ative magnification bias, b/bad = does not correspond 
exactly to the null line.) As the redshift distribution is 
peaked at z — 1.4, b/b^d can be interpreted as approx- 
imately parametrizing the amplitude of the bias at red- 
shift 1.4. From this interpretation we obtain a value of 
the bias at z w 1.4 of b ~ 2.5±0.6, which is in good agree- 
ment with previous measurements from quasar clustering 
(3lj . We can associate this bias with a host halo mass 
M 2 oo'- using the bias model of J32[, we obtain a halo mass 
of logm (M^on/Mpt) = 12.9^Q 5, consistent with previous 
estimates (To. l3l|. 

V. TESTING THE POWER SPECTRUM 

A. Null Tests 

We check our result and our pipeline with a number 
of null tests. In a simple first test, we cross-correlate the 
quasar distribution in one part of the sky with the lens- 
ing convergence in another; as seen in Fig. 5, the results 
are consistent with null as expected, with \ 2 — 6-5 for 
10 degrees of freedom for a fit to null. A more sophis- 
ticated test is to calculate the cross-correlation of the 
quasar maps with the curl component of the lensing de- 
flection (this differs from the convergence reconstructed 
earlier, which is gradient-like). The reconstructed curl 
map is expected to be zero (though it should contain re- 
construction noise), and hence the cross-correlation with 
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FIG. 4. Likelihood as a function of quasar bias divided by the 
fiducial bias, b/bad (we assume that the shape of the redshift 
dependence is constant and has the fiducial form of Fig. 3, 
and modify the amplitude of the bias function to calculate 
this likelihood.) Interpreting our measurement of b/bad = 
1.02 ± 0.24 as a bias at z m 1.4 (the peak in the quasar 
distribution), we obtain b = 2.5 ± 0.6 at this redshift. 

the quasar maps should be zero as well. We reconstruct 
the curl component of the estimator as in (though 
keeping the normalization and filters unchanged from the 
earlier convergence reconstruction, and simply replacing 
the dot product in /(L, 1) of Eq. 7 with a cross prod- 
uct projected onto the l x x l y direction). The cross- 
correlation of the lensing curl component with the quasar 
maps is also shown in Fig. 5. The error bars are calcu- 
lated from theory as before. As expected, this test is con- 
sistent with a null result, with a (somewhat low) value of 
X 2 = 3.5 for 10 degrees of freedom. 

B. Estimating Potential Systematic Contamination 

We estimate the magnitude of what are expected to be 
the largest contaminants: infrared (IR) sources, SZ clus- 
ters and galactic cirrus, which contribute flux to the CMB 
temperature maps. (The level of radio source power is 
much smaller as we can resolve and mask such sources 
down to low flux levels.) Contamination of the cross- 
power spectrum is conceivable because the sources of IR 
and SZ signal trace the underlying matter field as the 
quasars do, and galactic cirrus could reduce observed 
quasar counts by extinction. As explained in [8j, any 
IR or SZ contamination would appear as a negative bias, 
as a large IR or SZ signal in a certain direction increases 
local gradients; the lensing estimator falsely interprets 
this as a signature of demagnification of the CMB and 
hence estimates a spurious underdensity in this direction. 

To obtain an estimate of the contribution of IR and 
SZ contamination to the measured cross-power spectrum, 
we construct simulated quasar maps which are correlated 
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FIG. 5. Two successful null tests, both consistent with zero. 
Upper panel: the cross-power spectrum of quasar and lens- 
ing maps covering different parts of the sky (permutation null 
test). Lower panel: the cross-power spectrum of the recon- 
structed curl component of the lensing signal with the quasar 
maps (curl null test). 

with the IR and SZ maps of the simulations of [28[ . The 
quasar maps are constructed by randomly populating all 
halo positions (listed for all halos with masses greater 
than 6.8 x 10 12 M Q in the catalog supplied with (28j ) 
with quasars. In populating the halos we use a redshift- 
dependent probability of occupation such that the final 
simulated quasar map has the correct redshift distribu- 
tion. Cross-correlating this quasar map with the true 
convergence maps of the same simulation, the signal is 
consistent with the theory line of Fig. 2, which confirms 
that our simple simulation has approximately the correct 
bias (this is due to the average mass in the halo catalog 
being similar to the typical halo mass of a linear bias 
model consistent with the fiducial model). 

The level of systematic contamination in our estimator 
can be obtained from these simulated quasar, IR, and SZ 
maps. Keeping the same filtering and normalization as in 
our data, we reconstruct the IR contaminant to the lens- 



ing signal by replacing the temperature maps in Eq. 7 
with the simulated IR maps (which we rescaled using an 
appropriate factor as in [7] to match more recent con- 
straints on IR source flux). The cross-spectrum of the 
resulting map with the simulated quasar maps gives a 
negative spurious signal which is ~ 7% of the theoretical 
prediction for the lensing-quasar cross-correlation. Re- 
peating this analysis with the thermal SZ simulations 
gives a similar negative contamination of order w 5% of 
theory. (The analysis should overestimate the contam- 
ination, as the simulated quasars are placed exactly in 
the centers of the same halos that source the SZ and IR 
signal, neglecting any mis-centering effects.) Any sys- 
tematic error in our measurement of the lensing-quasar 
cross-power spectrum due to contamination from both IR 
and SZ sources should thus be significantly smaller than 
the size of the statistical error. In addition, the fact that 
these contamination signals are negative means that our 
detection of a positive lensing-quasar cross-power spec- 
trum cannot be due to such systematics. 

Finally, we estimate the level of contamination from 
galactic cirrus using the dust maps of (34[. We subtract 
a map of the signal at f48GHz induced by the dust (ob- 
tained from |28|) from the ACT temperature data, recon- 
struct lensing, and re-estimate the lensing-quasar cross- 
power spectrum. We find that the change in the cross- 
power spectrum is very small, of order 3% of the theo- 
retical prediction for the lensing-quasar cross-spectrum. 
(This is unsurprising, as in our analysis large scale power 
below I = 500 has been filtered out of the temperature 
maps.) Contamination of the cross-power by galactic cir- 
rus is thus negligible. 

VI. SUMMARY AND CONCLUSIONS 

In this work we measure the cross-correlation between 
ACT CMB lensing maps and maps of the quasar distri- 
bution made from the SDSS-XDQSO catalog. We de- 
tect the cross-power spectrum at 3.8cr significance, di- 
rectly confirming that quasars trace mass. We check our 
detection with null tests including a cross-correlation of 
quasars with the reconstructed curl component of lensing, 
which is found to be null as expected. Potential system- 
atic contamination is estimated and found to be negli- 
gible. From our detection we estimate the quasar bias. 
We measure b/bad = 1-02 ± 0.24; interpreting this as a 
bias at z ~ 1 .4 (the peak in the quasar distribution), we 
obtain b — 2.5 ±0.6 at this redshift (which corresponds to 
a host halo mass of log lo (M 2 oo/M ) = 12.9±oJ). Unlike 
measurements from clustering, this lensing measurement 
involves a direct comparison of the quasar distribution 
with the mass distribution, with little modeling required. 

The study of high-redshift mass tracers with CMB 
lensing is a new field. In the next few years, the signal- 
to-noise ratio on lensing measurements should improve 
by an order of magnitude with data from experiments 
such as Planck, ACTPol and SPTPol ACTPol in 
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particular should provide high signal-to-noise ratio lens- 
ing measurements which have considerable overlap with 
SDSS quasar fields. Higher signal-to-noise will allow con- 
straints on quasar biases as a function of redshift, lumi- 
nosity, color or other properties and will thus provide a 
wealth of information on the properties of quasars and 
the halos that host them. More precise bias measure- 
ments of both quasars and galaxies will also allow tests 
of dark energy properties [39[ and modified gravity [38[ . 
This work lies at the beginning of an exciting research 
program: the study of astrophysics and cosmology with 
CMB lensing cross-correlations. 
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